Boundary induced phase transitions in driven lattice gases with meta-stable states 
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We study the effect of met a- stability onto boundary induced phase transitions in a driven lattice 
gas. The phase diagram for open systems, parameterized by the input and output rates, consists of 
two regions corresponding to the free flow and jammed phase. Both have been entirely characterized. 
The microscopic states in the high density phase are shown to have an interesting striped structure, 
which undergoes a coarsening process, and survives in the thermodynamic limit. 



Driven lattice gases (DLG) are characterized by a non- 
vanishing mass-flow in the stationary state which is gen- 
erated by non-potential forces |lj . This feature of DLG 
has far reaching consequences, because even the station- 
ary state of the system is not described in the framework 
of standard equilibrium statistical mechanics J2|. In re- 
cent years the study of various non-equilibrium lattice 
gas models, which have important applications, e.g. in 
biological context H or as models for traffic flow 
led to a profound understanding of some key features of 
systems far from equilibrium. 

One of the most fascinating effects observed for this 
kind of models are so-called boundary induced phase 
transitions, i.e. the change of the bulk properties due to 
a variation of the boundary conditions . Such phase 
transitions have been extensively studied for a class of 
ID DLG with a continuous flow density relation for peri- 
odic boundary conditions (PBC). For some models even 
exact results for the stationary state exist [p|-|i"cfl . A phe- 
nomenological theory |ll| was recently able to predict the 
phase diagram in the case of open boundary conditions 
(OBC) also for more complicated models JI3]. However, 
less is known for DLG which exhibit meta-stable states for 
PBC [|l3|, although they are relevant in different fields. 
Here, we study one of the simplest DLG model which 
shows meta-stability for PBC. We first discuss briefly the 
case of PBC, but our main focus is to study the effect of 
meta-stability in the case of OBC. Using a combination of 
analytical and numerical methods we explore the struc- 
ture of the phase diagram and characterize the typical 
behavior in the different phases. 

We study a DLG model defined on a chain of length 
L. Each site of the chain may either be empty (x) or 
occupied with a particle of velocity zero (0) or one (I). In 
a first step all particles which moved in the previous time 
step keep their velocity if the next site is empty. With 
probability qo, velocity one is assigned to particles with 
velocity zero and an empty site ahead. In all other cases 
velocity zero is assigned to the particle. In the second 
step particles with velocity one move to the next site. 
Both steps are applied synchronously to all particles. 

This model may be interpreted as a special case of 
the recently introduced VDR model for traffic flow [pjj, 



i.e. we restrict the maximal velocity of the particles to 
one. For qo = I we recover the deterministic asymmetric 
exclusion process with parallel dynamics [fLof . 




FIG. 1. Flow density relation for PBC with qo = 0.25 and 
L — 1,000. Empty circles show simulation results for an 
initialization where all particles have velocity one and at least 
one empty site in front. Filled circles show results starting the 
simulation with a compact cluster of particles. These results 
are compared to our estimate for L — > oo (solid line). 

Following fllf] , the structure of the flow-density rela- 
tion, i.e. the fundamental diagram, in the case of PBC 
is a key to understand the behavior of the open chain. 
Fig. [j] shows the numerically established fundamental di- 
agram for our model, in comparison with: 

J{p) = p if < p < 0.5 

J(p) = q (l-p) if <? /(l + qo) < P < 1 

The form of the fundamental diagram can be easily un- 
derstood. For densities < p < 0.5 and suitable initial 
configurations (the initial velocity of the particles is one 
and no particle is blocked) all particles move determin- 
istically with velocity one. At larger densities p > 0.5 
the number of particles exceeds the number of holes, i.e. 
one cannot avoid the formation of clusters of particles. 
The simulation results indicate that in the steady state 
one observes phase separation, i.e. the formation of a 
large compact cluster of particles which coexists with 
a free-flow regime. The density in the free-flow regime 
is determined by the frequency of particles leaving the 
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cluster. The typical time the first particle needs to sep- 
arate from the cluster is simply given by T = q^ 1 ■ This 
time already determines the average gap between two 
particles in the low-density regime. Therefore the av- 
erage density, pf ree in the free-flow regime is given by 
Pfree = 1 / (T + 1) = qo / (1 + qo)- The number of mov- 
ing particles and therefore the flux follows from particle 
conservation (see also Q). 

Both solutions for J(p) coexist for pf ree < p < 1/2. 
For finite systems it is possible that the large cluster dis- 
solves due to fluctuations. In free flow states, fluctuations 
are absent and therefore no spontaneous cluster forma- 
tion takes place at densities p < 0.5. However, in the 
thermodynamic limit, any random initial configuration 
leads to a jammed stationnary state. 

After this short discussion of the system with PBC, we 
proceed with the open chain. OBC are implemented as 
follows : at the left boundary particles with velocity one 
are introduced with probability a if the first site of the 
chain is empty. At the right end the particles leave the 
chain from site L with probability (3, irrespective of their 
velocity. 

The values of a and (3 determine for given values of L 
and qo the bulk behavior of the system. First we notice 
that for (3 = 1 the only stochastic element is due to the 
left boundary, i.e. we recover the deterministic asymmet- 
ric exclusion process. Therefore the exact result for the 
flow JhJ is given by J(a) = a/(l + a). This result for the 
flow is expected to be valid not only for f3 = 1 but in the 
whole low density phase, because for low bulk densities 
the flow is controlled by the input of particles. 

Next, we consider the case of large bulk densities, i.e. 
we study the system's performance for large values of a 
and (3 « 1. For this case we expect to find frequently 
a compact cluster of particles at the right boundary. 
Therefore the mean time interval for a particle leaving the 
system is the sum of the average waiting time Tl = 1/(3 
of the last particle at site L and the average time which 
passes until the last site is occupied again, Tj = 1/qo- 
We find immediately that for large awl and (3 <C 1 the 
flow is given by 



1 = qoP 
Tl+Tj q +(3~ 



(1) 



This simple scenario is only valid if the density of holes 
fed into the system is very low. For larger values of (3 we 
expect that it is also possible that the cluster separates 
from the right boundary. In this case the density close to 
the exit is determined by the outflow from a cluster. The 
temporal headway t between two consecutive particles is 
distributed via Ph(t) = qo(l — <?o)* 1 where t > 1. Now 
we calculate the typical amount of time which is needed 
from the arrival of a particle (at t = 0) to the arrival of 
the second particle at (t = T) at site L. Here we have to 
distinguish two cases: (i) the particle arrives at site L — l 
before the first particle has left the system or (ii) that it 
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FIG. 2. Comparison of the estimated flow with simulation 
results. The simulation results have been obtained for a chain 
of length L = 500 and go = 0.25. (a) The Figure shows the 
perfect agreement of simulation results and the estimate J (a). 
We used /3 = 1. and (3 = 0.5 as output probabilities, (b) f3 
dependence of the flux (a = 0.5). The estimated flux and 
the simulation results agree for sufficiently small values of f3. 
Simulations for different lengths of the chain indicate that the 
deviations for (3 near one are due to finite size effects. 



arrives later, such that it reaches the last site without 
being blocked by the particle ahead. The probability 
that a particle gets blocked if it arrives at t is given by 
Pb(t) = (1 — /?)* and that it can pass without blockage 

byi 5 /(t) = EU/3( 1 -/3)*" 1 - 

This leads to the average value of T given by 



p f (t)(t + i) + p b (t)(t + ^ + -) 
P qo 



p qo 



(2) 



Also for this scenario we recover the same value (|l|) of 
the flow. These two kind of configurations are generic if 
the output probability controls the capacity of the sys- 
tem, i.e. in the presence of clusters. Therefore we expect 
to find a unique functional dependence of the flow, if the 
system is controlled by the capacity of the exit. The tran- 
sition between the a or f3 controlled parameter regime is 
expected to be located at 
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where both fluxes agree. 

Comparison of our prediction for the flow with simula- 
tion results shows excellent agreement (see fig. ||). Sim- 
ulations for different lengths of the chain show that the 
differences between our predictions and the simulations 
systematically vanish for larger system sizes. The origin 
of these finite size effects, i.e. the enhancement of the flow 
for large (3, is the possibility to find occasionally cluster- 
free micro-states. We indicated this by the broken lines 
in the phase diagram (Fig. ^). The absence of maximal 
current phase in the thermodynamic limit is related to 
the non-analyticity of the fundamental diagram [O . 




FIG. 3. Typical space-time plots in the striped phases. 
The figure shows the first ~ 400 sites of a chain of length 
L — 1, 000. (a) For SPI the size of the last clusters decreases 
rapidly close to the entrance, (a = 0.1, f3 = 0.1, go = 0.4) (b) 
For SPII the efficient particle reservoir at the entrance of the 
system augments the size of the leftmost cluster. 
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FIG. 4. Phase diagram of the open system for go = 0.25. 
The phases are in order: low density phase (LD), striped 
phase I (SPI), and striped phase II (SPII). At the line a — go 
one observes a constant density profile. The dashed lines show 
parameter regimes where in finite systems a systematically 
larger flux has been observed. 

We proceed characterizing the typical microscopic 
structure of the stationary state, in different parameter 



regimes. We distinguish three different domains, i.e. a 
low density domain (LD) for f3 > (3t(&) (LD), and two 
high density domains for (3 < (3t(a). Both high density 
domains correspond to a striped phase (SP), but they ex- 
hibit different behaviors in the vicinity of the entrance, 
depending whether a < qo (SPI) or a > qo (SPII). In the 
LD phase we observe the typical two domain structure 
known from the asymmetric exclusion process. Clusters 
which eventually separate from the high density domain 
at the right dissolve in a few time steps. This dissolution 
leads to the localization of the high density domain at the 
right. The system takes a bulk density pb{a) = J (a). 
At site L the density is given p^ = J(a)/f3. Fits of 
several density profiles follow the pure exponential form 
P(n) = Pb{a) + [pl - Pb(a)} exp(-(L - n)/g). 

In the high density region (SPI + SPII), clusters are 
able to reach entrance. This region can be described by 
considering that at the right boundary, some free-flow 
segments are injected within a compact cluster. Typi- 
cally, the free-flow segments reduce to isolated holes for 
(3 < qo, whereas for [3 > qo, they contain both empty 
sites and particles with velocity one. Both the free- 
flow segments and the compact clusters separating them 
move backwards, producing a striped structure in spatio- 
temporal plots (see fig. ||). As long as it is far enough 
from the left boundary, a cluster has the same in- and 
output rates, thus its width follows a non-biased random 
walk, whose average is a constant. However, there is a 
non-vanishing probability that the width of the cluster 
crosses the zero value, leading to the coalescence of the 
two neighboring free-flow segments, which will never be 
able to separate again afterwards. The excess particles 
are distributed among the neighboring clusters. Thus, 
as the stripes move towards the entrance, their widths 
increase. However, the average density remains a con- 
stant equal to the density selected on the left bound- 
ary p — J([3)/[3. When the clusters arrive near the left 
boundary and are exposed to the entrance flux, their 
width suddenly obeys a biased random walk. If a > qo 
(phase SPII), the width increases on average and the 
high input rate stabilizes the cluster until it reaches the 
boundary (fig. ||-b). If clusters and free-flow segments at 
a given position were monodisperse, then the sudden rate 
change would induce a linear density profile on the left, 
with a crossover to the constant bulk value occurring at 
site h = atf/(l + a) where tf is the average time in- 
terval during which free-flow is observed in the first site, 
i.e. the typical time between the disappearance of a jam 
and the arrival of the next. Due to the coarsening of 
the striped pattern in upstream direction one observes 
a subextensive growth of tf with increasing system size. 
The density in the first site is determined by noticing that 
the flux (0) is also equal to j = a(l — p\). If we measure 
tf numerically, we find a quite good agreement between 
our estimate and the measured density profile (see fig. |^) , 
except that the crossover region is in fact very large due 
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to size dispersion of clusters and free-flow segments. 
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FIG. 5. Density profile in the SPII region. The numerical 
result (solid line) is compared with the theoretical prediction 
(dashed line) based on an assumption of monodispersity and 
parameterized by tf which has to be determined numerically. 

If a < qo (phase SPI), the width of the leftmost cluster 
decreases on average due to the reduced capacity of the 
entrance. Some clusters will even not be able to reach the 
entrance, as seen on fig. [| As long as we are far from the 
transition line /3{(a), still, clusters are on average large 
enough so that most of them reach the entrance and the 
depletion region remains localized on the left. The cal- 
culation done in SPII for the density profile is still valid. 
The only difference is that now the linear part has a pos- 
itive slope. However, as we approach the transition line 
(3t (a) , more and more clusters do not reach the boundary. 
When a given cluster disappears before reaching the left 
end, the next one is exposed earlier to the low incoming 
flux, and thus the depletion region may invade the whole 
system. The line separating the depletion region from 
the constant density region becomes delocalized over the 
whole system as the transition line is approached. 

To conclude, we analyzed a DLG exhibiting meta- 
stable states. The structure of the phase diagram was 
explored by simulations as well as by a phenomenologi- 
cal description. The analytical predictions for each phase 
agree extremely well with the simulation results. The 
most characteristic feature is the spontaneous formation 
of simultaneously existing clusters for restricted outflow. 
The coarsening of these clusters as they flow backwards 
is due to the stochasticity of the model. We stress that 
this striped pattern survives in the thermodynamic limit. 
We expect that this behavior and more generally the 
structure of the phase diagram are generic for DLG with 
branched fundamental diagrams - this branched struc- 
ture being a signature of meta-stability. Moreover we 
established, in the limit of large system size, the absence 
of a maximal current phase. 

Besides a phenomenological similarity with the clog- 
ging phenomena in granular flow |]l6| , the practical rel- 
evance can be illustrated by the example of traffic flow, 
where the importance of boundary induced phase tran- 



sitions is well accepted |T4j . Although we chose a rather 
simplistic model and do not aim to give a complete de- 
scription of traffic flow, our results can be related to ex- 
perimental studies that revealed the existence of parallel 
existing clusters on highways |l5| ]. These experimental 
results now can be directly related to the existence of 
meta-stable states. Moreover, our result shows that jams 
can be observed far upstream even when boundary in- 
duced. Another interesting feature are the important 
finite size effects, which probably can be used for a sys- 
tematic flow optimization p7| . 
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